152 


MONTHLY  WEATHER  REVIEW 


Volume  129 


Comparison  between  Wavenumber  Truncation  and  Horizontal  Diffusion 

Methods  in  Spectral  Models 

Peter  C.  Chu,  Xiong-Shan  Chen,  and  Chenwu  Fan 

Department  of  Oceanography,  Naval  Postgraduate  School,  Monterey,  California 


(Manuscript  received  8  November  1999,  in  final  form  22  May  2000) 

ABSTRACT 

Commonly  used  horizontal  diffusion  and  wavenumber  truncation  dealiasing  methods  in  spectral  models  are 
verified  using  the  National  Center  for  Atmospheric  Research  Community  Climate  Model  version  3.  For  the  same 
horizontal  grid  resolution,  time  step,  physical  processes,  boundary  conditions,  and  initial  conditions,  the  simulated 
climate,  using  the  horizontal  diffusion  alone  model,  is  better  than  that  using  the  wavenumber  truncation  method. 
In  comparison  with  the  observed  climate  data,  the  global  root-mean-square  of  simulated  January  monthly  mean 
500-hPa  geopotential  using  the  horizontal  diffusion  alone  model  is  25%  less  than  that  using  the  wavenumber 
truncation  model.  However,  for  the  same  spectral  resolution,  the  wavenumber  truncation  model  (high  horizontal 
grid  resolution)  leads  to  more  accurate  solutions  than  the  horizontal  diffusion  model  (low  horizontal  grid  res¬ 
olution). 


1.  Introduction 

Silberman  (1954)  introduced  the  spectral  method  in 
meteorological  modeling,  solving  the  nondivergent  bar- 
otropic  vorticity  equation  using  an  interaction  coeffi¬ 
cient  method  to  calculate  the  nonlinear  terms.  Platzman 
(1960),  Bear  and  Platzman  (1961),  Kubota  et  al.  (1961), 
and  Ellsaesser  (1966)  also  discussed  the  integration  of 
the  vorticity  equation  using  the  same  method.  This 
method,  however,  was  not  practical  for  operational  use 
due  to  the  considerable  computational  power  needed  to 
calculate  the  nonlocal  sums  of  the  nonlinear  terms.  This 
problem  was  solved  by  Eliasen  et  al.  (1970),  Orszag 
(1970,  1971),  and  Machenhauer  and  Rasmussen  (1972) 
by  using  the  transform  method.  The  method  involves 
the  inverse  transform  of  model  variables  from  spectral 
space  to  physical  space,  the  computation  of  nonlinear 
advection  terms  at  grid  points,  and  the  transform  of  the 
nonlinear  terms  to  spectral  space. 

The  nonlinear  instability  was  first  studied  by  Phillips 
(1959),  who  found  that  the  cause  of  this  instability  is 
aliasing.  The  short  wavelength  of  the  waves  to  be  rep¬ 
resented  on  the  grid  is  2A  (where  2A  is  the  grid  size), 
but  the  waves  generated  by  the  nonlinear  terms  (from 
the  advection  term)  may  have  wavelengths  of  less  than 
2A.  Thus,  the  generated  waves  with  wavelengths  of  less 
than  2A  must  be  falsely  represented  (or  aliased)  by  the 
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waves  with  longer  wavelength  (of  which  most  are  2A 
or  4A).  This  computational  phenomenon  is  referred  to 
as  aliasing.  Repeated  aliasing  over  many  times  may  give 
rise  to  rapid  growth  of  wave  energy  in  wavelengths  2A 
and  4A  and  causes  the  nonlinear  instability. 

The  nonlinear  instability  occurs  not  only  in  spectral 
models  but  also  in  finite-difference  models  and  was 
originally  found  in  a  finite-difference  model  by  Phillips 
(1959).  Phillips  suggested  that  the  nonlinear  instability 
could  be  controlled  by  periodically  removing  wave¬ 
lengths  4A  and  smaller  by  using  a  Fourier  analysis,  or 
by  application  of  a  space  filter,  or  by  inclusion  of  an 
artificial  horizontal  diffusion  term.  Several  methods 
were  proposed  to  avoid  the  nonlinear  instability:  1)  the 
horizontal  diffusion  method  (Phillips  1959),  2)  the 
wavenumber  truncation  method  (Orszag  1970,  1971; 
Machenhauer  and  Rasmussen  1972),  and  3)  the  enstro- 
phy  conserved  difference  scheme  (Arakawa  1966).  The 
first  two  methods  are  commonly  used  in  spectral  models. 
The  horizontal  diffusion  method  uses  a  sufficiently  larg¬ 
er  value  to  suppress  the  nonlinear  instability.  The  wave- 
number  truncation  method  uses  the  3 M  +  1  point  dis¬ 
crete  Fourier  transform  in  the  transform  method  to  ob¬ 
tain  the  nonlinear  term  in  the  2 M  +  1  point  discrete 
Fourier  transform,  where  M  is  the  maximum  truncated 
wavenumber  of  the  solution.  In  fact,  this  method  is  to 
remove  or  to  truncate  the  high  wavenumbers  from  the 
solution  to  prevent  the  nonlinear  instability.  In  current 
spectral  and  finite-difference  models,  several  of  the 
above  methods  are  used  at  the  same  time,  or  at  least 
one  of  the  above  methods  is  used  to  suppress  the  non¬ 
linear  instability. 
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In  current  spectral  models  such  as  the  National  Center 
for  Atmospheric  Research  (NCAR)  Community  Climate 
Model  version  3  (CCM3)  and  the  European  Centre  for 
Medium-Range  Weather  Forecasts  (ECMWF)  model, 
both  the  horizontal  diffusion  method  and  the  wave- 
number  truncation  method  are  used  to  prevent  the  al¬ 
iasing  and  nonlinear  instability.  In  general,  if  the  wave- 
number  truncation  method  is  deleted  in  the  current 
CCM3,  the  larger  coefficients  of  the  horizontal  diffusion 
need  to  be  used  to  prevent  the  aliasing,  and,  on  the 
contrary,  if  the  horizontal  diffusion  method  is  deleted, 
a  stronger  wavenumber  truncation  method  needs  to  be 
used  in  the  model.  One  problem  emerges:  What  happens 
if  only  one  method  (horizontal  diffusion  or  wavenumber 
truncation)  is  used  in  the  spectral  model?  We  will  study 
this  problem  here  using  the  CCM3.  The  outline  of  the 
paper  is  as  follows.  In  section  2,  we  describe  the  climate 
spectral  model  using  CCM3.  In  sections  3  and  4,  we 
discuss  the  performance  of  the  wavenumber  truncation 
method  only  and  the  horizontal  diffusion  method  only, 
respectively.  In  section  5,  we  compare  between  the  two 
dealiased  methods  in  climate  simulation.  In  section  6, 
we  present  the  conclusions. 

2.  Climate  spectral  models 

a.  Model  description 

The  CCM3  evolved  from  the  Australian  spectral  mod¬ 
el  described  by  Bourke  et  al.  (1977)  and  McAvaney  et 
al.  (1978).  The  CCM3  uses  both  wavenumber  truncation 
and  horizontal  diffusion  methods  to  suppress  the  non¬ 
linear  instability.  The  wavenumber  truncation  method 
is  used  to  remove  the  highest  one-third  wave  compo¬ 
nents.  The  horizontal  diffusion  method  is  used  for  tem¬ 
perature,  vorticity,  and  divergence  with  the  default  val¬ 
ues:  e„  =  0.06  for  the  time  filter,  X(02)  =  0.25  X  106  m2 
s  1  for  the  horizontal  diffusion  in  a  linear  V2  form  in 
the  top  three  levels  of  the  model,  and  Kff  =  1.0  X  1016 
m4s~'  for  a  linear  biharmonic  V4  form  at  all  other  levels. 
The  time  step  in  this  study  is  taken  as  Ar  =  20  min. 

b.  Initial  and  boundary  conditions 

The  initial  data  used  in  this  study  are  for  the  1  Sep¬ 
tember  climatology  of  the  atmospheric  and  surface 
fields,  provided  by  the  NCAR  Climate  and  Global  Dy¬ 
namics  (CGD)  Division.  The  initial  data  include  tem¬ 
perature  T,  zonal  wind  component  u,  meridional  wind 
component  v,  water  vapor  specific  humidity  q,  surface 
geopotential  </>  ,,  surface  pressure  ps,  and  four  subsur¬ 
face  temperatures  Ts. 

The  time-invariant  boundary  data  are  surface  albedo, 
orography,  surface  aerodynamic  roughness  over  land, 
surface  evaporation  factor,  vegetation  type,  January 
snow  cover,  and  July  snow  cover.  The  time-varying 
boundary  data  include  the  sea  surface  temperature, 
which  contains  12  months  of  data  as  12  time  samples, 
and  the  ozone  volume  mixing  ratios. 


c.  Nomenclature 

We  use  TXXGYY(e,  K(T>,  Km)  to  identify  each  model 
run  with  the  maximum  wavenumber  (XX),  the  grid  points 
around  a  latitude  cycle  (YY),  the  time  filter  (e),  and 
horizontal  diffusion  coefficients  (Ka\  K,Ay).  For  example, 
the  CCM3  spectral  model  T21G64  with  default  values 
of  the  time  filter  (e)  and  horizontal  diffusion  coefficients 
(X(2),  Km)  can  be  represented  by  T21G64(e0,  X(02),  X(04)). 
It  has  32  Gaussian  latitudes  from  the  South  Pole  to  the 
North  Pole  and  18  vertical  levels  with  a  top  at  2.917  hPa, 
triangular  truncation  with  the  maximum  wavenumber  2 1 , 
and  64  grid  points  around  a  latitude  circle  (approximately 
a  5.6°  X  5.6°  transform  grid  and  64  X  32  X  18  grid 
points). 

d.  Experiment  design 

We  design  the  following  experiments  to  verify  the 
wavenumber  truncation  method  and  the  horizontal  dif¬ 
fusion  method  in  spectral  modeling:  1)  the  wavenumber 
truncation  method  only  (i.e.,  K(2)  =  Kw  =  e  =  0),  2) 
the  horizontal  diffusion  method  only  (keeping  all  the 
wave  components),  and  3)  both  the  wavenumber  trun¬ 
cation  and  horizontal  diffusion  methods  present  but  with 
different  spectral  resolutions. 

For  the  type-1  experiments,  we  integrate  the  two  spec¬ 
tral  models: 

W1  =  T21G64(0,  0,  0)  and  W2  =  T16G64(0,  0,  0). 

For  the  type-2  experiments,  we  integrate  the  three  spec¬ 
tral  models: 

HI  =  T31G64(e0,  Xf2*,  K §*>), 

H2  =  T31G64(e0,  2 Kg>,  2 X{,4>),  and 

H3  =  T31G64(e0,  7X<2>,  7X<4>). 

For  the  type-3  experiments,  we  integrate  the  three  spec¬ 
tral  models: 

R1  =  T21G64(e0,  X<2>,  X<4>), 

R2  =  T21G64(e0,  7X<2>,  7X<4>),  and 
R3  =  T31G96(e0,  X<2>,  X<4>). 

e.  Model  output 

1)  Global  statistics 

The  CCM3  model  computes  the  global  statistics  at 
each  time  step  for  diagnostic  purposes:  root-mean- 
square  (rms)  divergence,  RMSD  (s  ');  rms  absolute  vor¬ 
ticity,  RMSZ  (s_1);  rms  temperature,  RMST  (K);  mass 
integral,  STPS  (Pa);  moisture  integral,  STQ  (kg  m  2); 
the  maximum  Courant  number  for  the  horizontal  ve¬ 
locity  field,  HOR  (dimensionless);  the  maximum  Cour¬ 
ant  number  for  the  vertical  velocity  field,  VERT  (di¬ 
mensionless). 
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Fig.  1.  Temporally  varying  global  statistics  for  type-1  (wavenumber  truncation)  experiments  (a)  W1  and  (b)  W2:  RMSD  (in  10  5  s  *), 
RMSZ  (in  10~4  s_l),  RMST  (in  K),  STPS  (in  105  Pa),  STQ  (in  102  kg  m~2),  HOR  (dimensionless),  and  VERT  (dimensionless). 


2)  Monthly  mean  500-hPa  geopotential 

FIELDS 

Four  of  the  above  seven  spectral  models — H3,  Rl, 
R2,  R3 — are  selected  for  the  comparison  of  the  monthly 
mean  500-hPa  geopotential  fields.  We  integrate  the  four 
models  (H2,  Rl,  R2,  R3)  from  the  initial  climatological 
data  (1  Sep)  for  2  yr  and  compute  the  monthly  mean 
500-hPa  geopotential  fields  from  the  model  output. 

We  use  the  ECMWF  Reanalysis  sample  data  for 
1979-93  for  model  verification.  Since  the  four  models 
(H3,  Rl,  R2,  R3)  in  this  study  simulate  the  annual  var¬ 
iation  of  the  atmospheric  fields,  we  average  the  ECMWF 
Reanalysis  synoptic  sample  data  over  1979-93  to  obtain 
the  climatological  monthly  mean  data  for  verifying  the 
model  results. 

3.  Wavenumber  truncation  method 

To  verify  the  importance  of  the  wavenumber  trun¬ 
cation  method,  we  integrate  the  W1  and  W2  models  and 


show  the  results  in  Fig.  la  for  W1  and  Fig.  lb  for  W2. 
The  parameters  HOR  and  VERT  exceed  1.0  at  the  732d 
time  step  (equivalent  to  10  days)  in  the  W1  and  at  1023d 
time  step  (equivalent  to  14  days)  in  the  W2,  and  both 
models  blow  up.  This  indicates  that  use  of  the  wave- 
number  truncation  method  alone  cannot  suppress  all  the 
nonlinear  instability. 

4.  Horizontal  diffusion  method 

To  verify  the  importance  of  the  horizontal  diffusion 
method,  we  integrate  the  HI  and  H2  models  and  show 
the  results  in  Fig.  2a  for  HI  and  Fig.  2b  for  H2.  The 
parameters  HOR  and  VERT  in  the  HI  exceed  1.0  at  the 
144th  time  step  (equivalent  to  2  days),  and  the  model 
blows  up.  As  we  double  the  horizontal  diffusion  coef¬ 
ficients  (H2  model),  the  parameters  HOR  and  VERT  are 
always  less  than  1.0  (computationally  stable);  however, 
the  model  output  contains  high-frequency  noises. 
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Fig.  2.  Temporally  varying  global  statistics  for  type-2  (horizontal  diffusion)  experiments  (a)  HI  and  (b)  H2:  RMSD  (in  10-5  s_l),  RMSZ 
(in  10-4  s_1),  RMST  (in  K),  STPS  (in  105  Pa),  STQ  (in  101 2  kg  m-2),  HOR  (dimensionless),  and  VERT  (dimensionless). 


5.  Comparison  among  four  stable  models 

Following  is  the  comparison  among  the  four  stable 
model  runs:  H2,  Rl,  R2,  and  R3. 

a.  Monthly  mean  500-hPa  geopotential  fields 

Each  model  has  two  simulated  January  and  two  sim¬ 
ulated  July  fields.  We  compare  the  second  simulated 
January  and  July  with  the  monthly  mean  ECMWF  data 
(taken  as  the  “observed”  data).  The  comparison  is  con¬ 
ducted  in  two  subgroups:  (a)  same  grid  resolution  (H3, 
Rl,  R2),  and  (b)  same  spectral  resolution  (H3,  R3). 

1)  Models  with  the  same  grid  resolution 
(H3,  Rl,  R2) 

Figure  3  shows  the  observed  and  the  model-simulated 
monthly  mean  January  500-hPa  geopotentials.  The 
troughs  of  east  Asia  and  the  east  coast  of  North  America 
are  well  simulated  in  strength  and  in  location  in  the  H3 


model.  The  simulated  troughs  in  the  Rl  and  R2  models 
are  weaker  in  strength.  The  trough  of  the  Middle  East 
is  well  simulated  in  all  the  climate  models  in  comparison 
with  the  observed  data.  In  the  Southern  Hemisphere, 
the  straight  jet  stream  is  well  simulated  in  the  H3  model, 
but  the  jet  streams  in  the  Rl  and  R2  models  are  too 
weak  and  meandering. 

Figure  4  shows  the  observed  and  the  model-sim¬ 
ulated  monthly  mean  July  500-hPa  geopotentials.  In 
the  Northern  Hemisphere,  the  ridge  of  North  America 
is  well  simulated  in  the  three  models.  In  the  Southern 
Hemisphere  the  jet  stream  is  too  strong  in  the  H3 
model  and  is  too  weak  in  the  Rl  and  R2  models.  The 
observed  geopotential  field  shows  a  deep  low  with  a 
geopotential  contour  of  47  600  m2  s~2  near  the  South 
Pole.  The  H3  model  (horizontal  diffusion  alone)  sim¬ 
ulates  this  deep  low  very  well  with  the  same  contour 
(47  600  m2  s~2)  at  the  same  location  (near  the  South 
Pole).  However,  the  Rl  and  R2  models  fail  to  simulate 
this  system. 


156 


MONTHLY  WEATHER  REVIEW 


Volume  129 


Longitude 


Longitude 


Fig.  3.  The  observed  and  simulated  by  stable  models  (H3,  Rl,  R2, 
R3)  Jan  monthly  mean  500-hPa  geopotentials  (in  10  m2  s-2). 


Fig.  4.  The  observed  and  simulated  by  stable  models  (H3,  Rl,  R2, 
R3)  Jul  monthly  mean  500-hPa  geopotentials  (in  10  m2  s-2). 


2)  Models  with  the  same  spectral  resolution 

(H3,  R3) 

Both  H3  and  R3  in  January  have  the  same  spectral 
resolution  with  the  maximum  wavenumber  3 1 .  The  sim¬ 


ulated  troughs  of  east  Asia  and  the  east  coast  of  North 
America  in  the  R3  model  are  weaker  in  strength  than 
in  the  H3  model.  The  trough  of  the  Middle  East  is  well 
simulated  in  the  H3  model  but  not  in  the  R3  model.  In 
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Geopotentiais  Geopotentials 


Fig.  5.  The  Jan  and  Jul  meridional  cross  sections  of  the  observed  and  the  simulated  zonal- 
averaged  monthly  mean  500-hPa  geopotentials  (in  10  m2  s-2). 


the  Southern  Hemisphere,  the  straight  jet  stream  is  well 
simulated  in  both  the  H3  and  R3  models  with  compa¬ 
rable  strength  to  the  observed  field  (Fig.  3). 

In  the  Northern  Hemisphere,  the  ridge  of  North  Amer¬ 
ica  in  July  (Fig.  4)  is  well  simulated  in  both  models.  In 
the  Southern  Hemisphere,  the  jet  stream  is  too  strong 
in  the  H3  model  and  is  quite  reasonable  in  the  R3  model. 
The  observed  geopotential  field  shows  a  deep  low  with 
a  geopotential  contour  of  47  600  m2  s  2  near  the  South 
Pole.  The  H3  model  simulates  this  deep  low  very  well 
with  the  same  contour  (47  600  m2  s  2)  at  the  same  lo¬ 
cation  (near  the  South  Pole).  However,  the  R3  model 
fails  to  simulate  this  system. 

b.  Zonally  averaged  monthly  mean  500-hPa 

geopotentials 

We  compute  the  zonally  averaged  monthly  mean  500- 
hPa  geopotentials.  From  the  meridional  cross  sections 
of  the  zonal-averaged  monthly  mean  500-hPa  geopo¬ 
tentials  (Fig.  5),  we  find  that  among  the  same  grid  res¬ 
olution  models  (H3,  Rl,  R2)  the  H3  model-simulated 
values  are  much  closer  to  the  observed  data  than  the 
other  two  (Rl,  R2)  model-simulated  values,  especially 
in  the  middle  and  high  latitudes. 


Table  1.  Root-mean-square  error  of  Jan  mean  500-hPa 
geopotentials  (in  10  m2  s~2). 


Model 

Global 

Northern 

Hemisphere 

Southern 

Hemisphere 

H3 

216 

234 

199 

Rl 

289 

276 

300 

R2 

281 

275 

286 

R3 

106 

109 

104 

For  the  same  grid  resolution,  comparable  values  sim¬ 
ulated  using  Rl  to  using  R2  (same  wavenumber  trun¬ 
cation  but  different  horizontal  diffusion  coefficients)  im¬ 
plies  the  insensitivity  of  model  results  to  horizontal  dif¬ 
fusion  coefficients.  Improvement  of  the  simulation  (es¬ 
pecially  in  high  latitudes)  using  H3  compared  to  using 
R2  (same  large  horizontal  diffusion  coefficients  but  dif¬ 
ferent  spectral  resolution)  indicates  that  a  high  spectral 
resolution  (untruncated  model)  leads  to  a  high  model 
accuracy.  For  the  same  spectral  resolution,  the  high  grid- 
resolution  model  (R3)  simulates  values  closer  to  the 
observed  values  than  the  H3  model  although  it  con¬ 
sumes  much  more  computer  time.  Thus,  both  the  grid 
resolution  and  the  spectral  resolution  are  key  issues  in 
determining  the  spectral  model  accuracy. 

Tables  1  and  2  indicate  the  rms  errors  of  the  monthly 
mean  500-hPa  geopotentials  between  observed  and  sim¬ 
ulated  values  for  the  globe,  the  Northern  Hemisphere, 
and  the  Southern  Hemisphere,  respectively.  Among  the 
same  grid  resolution  models  (H3,  Rl,  R2),  the  horizontal 
diffusion  alone  model  (H3)  has  the  least  errors  in  sim¬ 
ulating  climatological  500-hPa  fields.  Such  an  improve¬ 
ment  is  more  evident  in  January  than  in  July.  In  January, 
the  rms  errors  of  the  H3  model  for  the  globe,  the  North¬ 
ern  Hemisphere,  and  the  Southern  Hemisphere  are  2160, 


Table  2.  Root-mean-square  error  of  Jul  mean  500-hPa 
geopotentials  (in  10  m2  s-2). 


Model 

Global 

Northern 

Hemisphere 

Southern 

Hemisphere 

H3 

228 

253 

203 

Rl 

260 

255 

265 

R2 

238 

258 

216 

R3 

103 

118 

87 

158 
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2340,  and  1990  m2  s~2,  respectively.  These  values  are 
a  25%,  15%,  and  34%  error  reduction  compared  to  that 
of  the  R1  model  and  are  a  23%,  15%,  and  30%  error 
reduction  compared  to  that  of  the  R2  model.  In  July, 
the  rms  errors  of  the  H3  model  for  the  globe,  the  North¬ 
ern  Hemisphere,  and  the  Southern  Hemisphere  are  2280, 
2530,  and  2030  m2  s~2.  These  values  are  a  12%,  1%, 
and  23%  error  reduction  compared  to  that  of  the  R1 
model,  and  are  a  4%,  2%,  and  6%  error  reduction  com¬ 
pared  to  that  of  the  R2  model,  respectively.  This  indi¬ 
cates  an  encouraging  improvement  in  the  climate  sim¬ 
ulation  using  the  horizontal  diffusion  alone  method 
among  the  same  grid  resolution  model.  For  the  same 
spectral  resolution  models  (H3  and  R3),  the  high  grid 
resolution  model  (R3)  leads  to  a  more  accurate  solution. 

6.  Conclusions 

Our  conclusions  are  as  follows. 

1)  The  capability  of  horizontal  diffusion  and  wave- 
number  truncation  dealiasing  methods  to  suppress 
nonlinear  instability  in  spectral  models  is  investi¬ 
gated  using  the  NCAR  CCM3  model.  The  wave- 
number  truncation  method  alone  cannot  suppress  the 
nonlinear  instability  even  half  of  the  wavenumbers 
are  truncated.  The  horizontal  diffusion  method  alone 
can  suppress  the  nonlinear  instability  when  the  hor¬ 
izontal  diffusion  coefficients  exceed  some  critical 
values. 

2)  When  the  horizontal  diffusion  coefficients  (W(2),  W(4)) 
are  taken  at  twice  the  CCM3’s  default  values  (0.25 
X  106  m2  s_1,  1.0  X  1016  m4  s'1),  the  horizontal 
diffusion  alone  model  is  computationally  stable,  but 
the  solution  contains  high-frequency  noises.  Increase 
of  the  horizontal  diffusion  coefficients  leads  to  fil¬ 
tration  of  high-frequency  noises. 

3)  Both  the  grid  resolution  and  the  spectral  resolution 
are  key  issues  in  determining  the  spectral  model  ac¬ 
curacy.  For  the  same  grid  resolution,  the  simulated 
climatological  monthly  mean  fields  are  relatively  in¬ 
sensitive  to  horizontal  diffusion  coefficients  and  are 
improved  using  an  untruncated  model  (i.e.,  highest 
spectral  resolution),  when  the  diffusion  coefficients 
are  sufficiently  large.  For  the  same  spectral  resolu¬ 
tion,  the  high  grid  resolution  model  leads  to  a  more 


accurate  simulation.  This  implies  the  importance  of 
the  wavenumber  truncation  method. 
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